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While the photon forward-scattering amplitude on free magnetic dipoles (e.g. free neutrons) van- 
ishes, the nucleon magnetic moments still contribute significantly to the photon dispersion relation 
in a supernova (SN) core where the nucleon spins are not free due to their interaction. We study the 
frequency dependence of the relevant spin susceptibility in a toy model with only neutrons which 
interact by one-pion exchange. Our approach amounts to calculating the photon absorption rate 
from the inverse bremsstrahlung process ynn — > nn, and then deriving the refractive index n rc f r 
with the help of the Kramers-Kronig relation. In the static limit (lu — > 0) the dispersion relation is 
governed by the Pauli susceptibility XPauii so that n 2 cfl — 1 w ^Pauii > 0. For lu somewhat above the 
neutron spin-relaxation rate r CT we find n 2 cfl — 1 < 0, and for u»r„ the photon dispersion relation 
acquires the form lu 2 — k 2 = m 2 . An exact expression for the "transverse photon mass" m 7 is given 
in terms of the /-sum of the neutron spin autocorrelation function; an estimate is rn 2 ~ XPauiiTTo-- 
The dominant contribution to n rc f r in a SN core remains the electron plasma frequency so that the 
Cherenkov processes 7^ <-> v remain forbidden for all photon frequencies. 

PACS numbers: 97.60.Bw, 95.30. Cq, 52.25.Tx, 14.60.Lm 



I. INTRODUCTION 



The photon dispersion relation in astrophysical plas- 
mas is usually dominated by the electromagnetic inter- 
action with the electrons of the medium. It was recently 
claimed [Q, however, that in a supernova (SN) core the 
photon interaction with the magnetic moments of the 
nucleons yields the dominant contribution to the refrac- 
tive index n re f r . Because the new contribution has the 
opposite sign of the usual plasma term, the photon four- 
momentum K would actually become space-like, allowing 
for the Cherenkov processes jv — ■* v and v — > 1/7 which 
could be of great importance for the neutrino opacities. 
Due to a numerical error in Ref. [jjj the overall magni- 
tude of the nucleon magnetic-moment effect is in fact 
much smaller 0, but even after this correction it is not 
very much smaller than the electron plasma effect and 
thus deserves a closer look. 

It is surprising, at first, that the nucleon magnetic mo- 
ments contribute at all to the refractive index because the 
photon forward-scattering amplitude on fermions with a 
magnetic moment is identically zero. Most recently the 
photon polarization tensor for an ensemble of noninter- 
acting spin-i particles was calculated in great detail ||] 
and it was found that indeed the magnetic moments alone 
do not produce any contribution to n re f r . However, the 
underlying assumption of a collisionlcss system is far from 



justified in a SN core where the nucleon spin-spin inter- 
action plays a dominant role. For photon frequencies 
below the nucleon spin relaxation rate T a the hydrody- 
namic limit is the physically appropriate description (not 
the collisionless limit) justifying the use of the Pauli sus- 
ceptibility in Ref. 0. In a SN core the spin relaxation 
rate is likely to be of the same order as the temperature 
T 0, typical photon frequencies are also of that order, 
so that for the entire spectrum of relevant photon fre- 
quencies neither the hydrodynamic nor the collisionless 
limits are truly justified. Therefore, an understanding of 
the photon refractive index and its frequency dependence 
requires a more general analysis than has been offered in 
either Ref. § or §. 

Perhaps the easiest way to appreciate this point is to 
consider photon absorption. In a collisionless system 
of neutral fermions ("neutrons") with a magnetic mo- 
ment \x n the refractive index is n rcfr = 1 up to order 
/i^ which implies that there is no Landau damping, i.e. 
no Cherenkov effect 771 <-* n. The only photon damp- 
ing occurs at order /i„ from magnetic Compton scatter- 
ing 771 — > 717. On the other hand, if our "neutrons" 
interact by a spin-dependent force (which for real neu- 
trons is caused by pion exchange) we have the inverse- 
bremsstrahlung absorption process 77m — > nn so that we 
do have photon absorption to order /i^. Its rate far ex- 
ceeds that of magnetic Compton scattering because in a 
SN core the nn interaction rate is large. By virtue of the 
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Kramers-Kronig relation one can then derive the associ- 
ated refractive index n re f r which does not vanish to order 
//„. (Note that we always take n re f r to be a real quantity 
even though one sometimes describes absorption by an 
"imaginary part of the refractive index.") 

We proceed in Sec. [El] with the general photon disper- 
sion relation in a pure neutron medium in terms of the dy- 
namical spin-density structure function by virtue of the 
fluctuation-dissipation theorem and the Kramers-Kronig 
relation. In Sec. Ill we use a semi-heuristic expression 
for the dynamical spin-density structure function in the 
long-wavelength limit to obtain a quantitative estimate 
of the magnetic-moment refractive index in a SN core. 
In Sec. IV we discuss and summarize our findings. 



II. GENERAL RELATIONS 



caused by an applied magnetic induction. However, the 
magnetic field H and the transverse part of the displace- 
ment, characterized by k • Dy = 0, do not have inde- 
pendent meaning Therefore, among other possibili- 
ties one may equally well choose H = B, = etEt, 
and Di = tj^L with ei, = e the longitudinal and 
€t = e+ (1 — k 2 /u} 2 the transverse dielectric permit- 
tivity. They are related to the polarization functions by 
(l = 1 - ttl/(w 2 - k 2 ) and e T = 1 - tt t /uj 2 @. In this 
language the dispersion relations take on their standard 
form e L (uj,k) — and w 2 e T (w, k) = k 2 0. 

We are presently only concerned with the disper- 
sion relation of transverse modes ("photons") because 
a medium consisting of magnetic dipoles is not expected 
to support longitudinal modes (longitudinal plasmons). 
From the above it follows immediately that the photon 
dispersion relation can be written in the form 



A. Photon Dispersion 

The main idea behind our treatment of the photon dis- 
persion in a neutron medium is the observation that the 
photon absorption rate r a b s is dominated by the inverse- 
bremsstrahlung process 77m — > nn which is enabled by 
the tensor component of the pion-exchange force between 
neutrons. From the absorption rate we can determine the 
refractive index n re f r with the help of a Kramers-Kronig 
relation. The inverse-bremsstrahlung process itself can 
be calculated easily by the usual perturbative methods. 
However, a SN core is so dense and hot that these meth- 
ods are not obviously justified |Q. Therefore, it is more 
useful to begin with a discussion of photon absorption in 
the language of linear-response theory which allows us to 
identify more general properties of the photon refractive 
index than would be apparent on the perturbative level. 
In order to apply our general results to a SN core we will 
then, of course, have to take recourse to a semi-heuristic 



perturbative approach to nn interactions (Sec. III). 

Photon dispersion is caused by the medium's response 
to applied electromagnetic fields. In the homogeneous 
and stationary case all relevant information is contained 
in the polarization tensor Xl^ v (K) where K — (w,k) is 
the frequency and wavevector of the applied electromag- 
netic perturbation. If the medium is isotropic and parity 
conserving the polarization tensor is uniquely character- 
ized by a pair of two response functions which are often 
chosen to be the longitudinal and transverse polarization 
functions n L = (1 - uj 2 /k 2 )U 0Q and n T = ±(Trn - tt l ) 
with k = |k| the wavenumber of the perturbation. Be- 
cause of the assumed isotropy all quantities depend only 
on the wavenumber k, not on the direction of k. The 
dispersion relation of propagating modes is determined 
by uj 2 — k 2 = 7Tt,l(w, k). 

Another pair of medium response function is the di- 
electric permittivity e and the magnetic permeability fi 
which give us the displacement D = eE caused by an 
applied electric field and the magnetic field H = /i _1 B 



(1) 



With the usual definition of the photon refractive index 

k 



7Z rc f r 



(2) 



we arrive at the classical result n 2 efr = e/i |. It fol- 
lows that the refractive index must be determined self- 
consistently as a solution of 



"?efr(w, k) = e(w, k) u(uj, k) 



(3) 



with k = n te f T uj for any frequency w of a propagat- 
ing mode. Depending on the properties of the medium 
the long- wavelength approximation e(aj,k) w e(w,0) and 
li(uj,k) ~ /i(w,0) may be justified, leading to the much 
simpler dispersion relation n 2 efr (ijj) = e(u>, 0) /u(k>, 0). 

Sometimes it will be more useful to write the photon 
dispersion relation in the form uj 2 — k 2 - 
of a frequency-dependent "effective mass 



TO^ff m terms 



'oil' 



H = (l-n 2 cfr )c 2 



(4) 



where in fact m 2 ff < is possible. For electric interac- 
tions and frequencies well above all resonances we obtain 
the well-known plasma effect dispersion relation which 
implies that m 2 ff > and independent of frequency [||. 
We will show that the same holds true for our magnetic 
case. Therefore, it is useful to define 



i 7 = lim m e fi(u>) 



(5) 



as the (transverse) "photon mass" in the medium. 

We will mostly be concerned with a medium of neu- 
trons which interact with the electromagnetic field by 
their magnetic dipole moment. In the nonrelativistic 
limit they do not respond at all to an applied electric 
field so that we may use the approximation e = 1. The 
magnetic permeability can be written as fi = 1 + x m 
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terms of the magnetic susceptibility \- (We use rational- 
ized units where a = e 2 /47r w 1/137 or else we would 
have to write fj, = 1 + 47rx ||.) 

In general the magnetic susceptibility is a complex 
function of the real variables uj and k. Following com- 
mon practice we write it in the form 



x(w,A)=xV>*) + *X // («',*) 



(6) 



In our case the most useful quantity to describe the 
neutron spin fluctuations is the dynamical spin-density 
structure function. Following the normalization conven- 
tion of Ref. it is defined as 



S a (uj, k) 



4 
3n r , 



dte iuJt ((r(t,-k)-(T(0,--k)), (9) 



in terms of its real and imaginary parts. It is found that 
that x" is an °dd function of uj while \' is even Q . Be- 
cause we have defined the refractive index to be a real 
quantity the dispersion relation is 



r (w,fc) - 1 = x'(w>*0 



(7) 



with k — n Te f T LO. 

The imaginary part of the susceptibility describes dis- 
sipation: Usually one pictures a stationary beam of fre- 
quency uj along the z-direction which is characterized by 
a (real) wavenumber k = n re f r w and a damping wavcnum- 



ber 



with A the photon mean free path. The 



amplitude of this beam varies as e 
sity as e 



-i(u)t— kz)- 



its inten- 



-2kz 



implies n^ eb 
1 



e z l x . The relativistic limit |n r cfr— 1| <C 1 

1 = ("reft- + l)(«-rcfr - 1) » 2(n re fr ~ 1) Or 

^rcfr — i ~ hx- Therefore, one can picture \\" to be an 
"imaginary part of the refractive index" , leading to the 
identification k = 5%"^ or x" = (Aw) -1 . 

We stress that at finite temperature this simple inter- 
pretation is not complete because the medium can both 
absorb and spontaneously emit photons. The two pro- 
cesses are related by the usual detailed-balance factor 
e~ u / T . What is actually damped is not a mode k of the 
electromagnetic field, but rather the deviation of its oc- 
cupation number from a thermal distribution. It is easy 
to show that this damping occurs at a rate 1 — e~ u ^ T 
times the "naive" absorption rate r a b s flC|] . Therefore, 
the appropriate interpretation of the imaginary part of 
the susceptibility is 



X "(u>,k) 



l-e 



J ' T ) r abs M 



(8) 



with k = n Te f T uj. In the limit |?i ro f r — 1| <C 1 the "naive" 
absorption rate is r a b s = A -1 ; it is given by the standard 
formula "absorption cross section times target density." 



B. Fluctuation-Dissipation Theorem 



where cr{t 1 k) is the spatial Fourier transform of the neu- 
tron spin-density operator er(t,r) = ^^(t, r)Tt(j(t, r). 
Here ip(t, r) is a Pauli two-spinor describing the nucleon 
field and r is a vector of Pauli matrices. Further, n n is 
the neutron number density and (. . .) denotes a thermal 
ensemble average. Of course, in an isotropic system the 
structure function depends only on k = |k|. The normal- 
ization was chosen such that 



(10) 



for a case where there are no static spin-spin correlations 
between different neutrons which are taken to be nondc- 
generate. In the limit of vanishing spin-spin interactions 
we have 



S a (oj,0) -> 2tt8{lo). 
Moreover, it satisfies 

S a (-oj,-k) = e- u,T S a (u),V) 



(11) 



(12) 



as required by the principle of detailed balance. 

We next observe that the operator for the magnetiza- 
tion density for neutrons is M — 2/x ra cr where the factor 
2 is the gyromagnetic ratio for a spin-i particle and fi n is 
the neutron magnetic moment, not to be confused with 
the magnetic permeability fi of the previous section. A 
relationship between the dissipative part of x an d sponta- 
neous spin fluctuations can then be written in the form 



X>,k) = - 



-foe 



dte iut 



1^ 



3 ^ 



[M^kJ.M^O.-k)] Y 
(13) 



To lowest order the neutrons can absorb photons only 
because they interact by a spin dependent force which en- 
ables the inverse-bremsstrahlung process jnn — > nn. At 
the same time this spin-dependent force causes the neu- 
tron spins to fluctuate. The relationship between spin 
fluctuations and the absorptive part of the spin suscepti- 
bility is encapsulcd in the fluctuation-dissipation theorem 
which will help us to understand some general properties 
of the photon refractive index. 



where [ • , • ] is the usual commutator. Comparing this 
with Eqs. (j^) and jl2] ) reveals that this relationship is 
equivalent to 

X >,k) = ~nln n (l - e-/ T ) S a (uj,k). (14) 



In this form it is known as the fluctuation-dissipation 
theorem MH 
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C. Kramers-Kronig Relation 

Once the imaginary part of the magnetic susceptibility 
is known the real part can be found by virtue of the well- 
known Kramers-Kronig relation 



, du> x"(&>k) 
X {uJ,k)=P —■ 



7T UJ — UJ 



(15) 



where P denotes a Cauchy principal value integral. With 
the help of the fluctuation-dissipation theorem Eq. ( |l4| ) 
we find a direct relationship between the dispersive part 
of the magnetic susceptibility and the spin-density struc- 
ture function 

'i u\ o 2 f +oc duj uj S a (uj,k) 
X {oj, k) = 2n n n n P / — — . (16) 

J~oo 27T UJ 2 - UJ 2 

One may use detailed balance to write this in the form 
'( ^ D [°°du l~e*/ T uJ 2 S a (uj,k) 

X (W, k) = XPauli P / T-f^ ^ o - , (17) 



q 7T Lj/T UJ 2 — LU 2 



where 



XPauli 



2 

T 



(18) 



is the usual Pauli susceptibility for a system of collision- 
less spin-i particles with a magnetic moment /i n . 



D. Limiting Cases 

In order to understand the general behavior of the re- 
fractive index we begin with the static limit uj — > 0. The 
static susceptibility xo(^) = x'(Oj^) nas 110 imaginary 
part because x" is an °dd function of uj. From Eq. ( |17| ) 
we find for the real part 



Xo(k) = XPauli 



duj 1 



o/T 



Cj/T 



S a (tu,k). (19) 



In the collisionless limit the structure function becomes 
narrowly peaked around uj — 0. With the help of Eq. ( fTT| ) 
we thus recover the usual long-wavelength result result 
Xo(0) = XPauli- When S a (uj,k) is not narrowly peaked 
on scales of the temperature, the static susceptibility de- 
creases relative to the Pauli value — we show this effect 
explicitly in Fig. ^ below in the framework of a heuristic 
toy model. 

How large may the frequencies be that the static result 
is still approximately justified? The structure function in 
the long-wavelength limit S a (uj, 0) has the interpretation 
of the autocorrelation function of a single neutron spin. 
Therefore, it is a broad, decreasing function of u> with 
a width representing something like the spin-relaxation 
or spin- fluctuation rate T a . If the external electromag- 
netic perturbation has a frequency much less than this, 



w < r 5 , we are in the hydrodynamic limit where the 
neutron spins may fully relax to a new thermodynamic 
equilibrium state on the time scale of a period of the per- 
turbation. In this case we may use the static suscepti- 
bility to estimate the photon refractive index. Moreover, 
even though we just saw that the static susceptibility is 
not independent of the width of S a (uj), this dependence 
is weak so that in the hydrodynamic limit the Pauli sus- 
ceptibility is a good estimate, justifying the approach of 
Ref. [Q to photon dispersion in the limit uj <C T^. 

The opposite limiting case is that of very large u). If 
S a {uj) falls off sufficiently fast beyond some frequency lu 
which is determined by the nature of the nn interaction 
potential, then for uj ^> luq the integral in Eq. ( |l6| ) is 
dominated by \Q\ ^ ujq, leading to 



X '(uJ,k) = -2 



^CjS a (^,k). (20) 



The integral in this equation is the so-called /-sum of 
the structure function. Independently of the nature of 
the assumed nn interaction the /-sum always exists and 
is given as a thermal expectation value of the tensor part 
of the nn interaction potential ]l^| . Moreover, the /- 
sum is always positive because of the detailed-balance 
property Eq. ( |l2| ) . 

For photon dispersion, this result corresponds to a 
positive value for the squared effective mass defined in 
Eq. §). With Eq. (0) we find 



— ujSo-(uj, k). 



(21) 



If the momentum dependence of this expression is weak 
so that we may use the long-wavelength limit then the 
photon dispersion relation is that of a massive particle 
uj 2 — k 2 = with the transverse photon mass given 
by Eq. ( ^l|) with k — on the right-hand side. The 
appearance of this form has the same cause as in the 
usual plasma case, i.e. n rc f r is given by the /-sum of the 
relevant dynamical structure functions. 

The Pauli susceptibility is a positive number (the neu- 
trons are a paramagnetic medium) so that in the hy- 
drodynamic limit the photon dispersion relation is ap- 
proximately characterized by n^ ch . — 1 = XPauli or m^ s = 
-XPauli^ 2 < 0. On the other hand in the large- frequency 

Moreover, 
<e on the ap- 



limit we have m 2 ff > as given in Eq 
on dimensional grounds the /-sum must ta 
proximate value r^. Therefore 



'off 



XPauli X 



— UJ 

+ TT„ 



for uj < r o 
for uj > r o 



(22) 



gives us a rough picture of the behavior of the photon 
dispersion relation in a medium of neutron spins. 
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III. SEMI-HEURISTIC MODEL 

In a SN core neither the collisionless nor the hydro- 
dynamic limits are appropriate so that we need to come 
up with a concrete expression for the dynamical spin- 
density structure function in order to estimate the pho- 
ton refractive index. In a dilute medium one may use 
the usual perturbative methods to compute the processes 
jnn «-» nn. Because the relevant photon energies are 
small compared with the neutron mass the momentum 
transfer of the radiation to the neutron system may be 
neglected, an approximation which amounts to the long- 
wavelength limit which we shall henceforth adopt with 
the notation Scr(u)) = lim/c_»o S a {pj, k). Next, one may 
extract (u>), where the superscript indicates that this 
is a lowest-order perturbative result. Independently of 
the details of the assumed nn interaction potential one 
finds the generic representation [Q 



S«(a,) = -| S (aVT) 



(23) 



where s(x) with x = uj/T is a slowly varying function of 
order unity. This factorization is somewhat arbitrary; we 
define what we call the neutron "spin-fluctuation rate" 
T a such that for nondegenerate neutrons s(0) = 1. More- 
over, we have 



s(—x) — s(x) e 



(24) 



so that the detailed-balance relation for S a (uj) Eq. fli2| ) 
is satisfied. 

The lowest-order perturbative representation sj^\uj) 
diverges at u> = and thus violates the normalization rule 
Eq. (|o|). However, including multiple-scattering effects 
suggests the "resummed" representation [Q 



r„ 



T 2 /4 



s(u/T). 



(25) 



In a very dilute medium this function is strongly peaked 
around lj = so that it approaches 2ttS(lo). In this limit 
we have r = T a , i.e. we approach the classical limit of a 
Lorentzian correlation function S a (u)) = Y a j (u 2 + T 2 /A). 
We stress that the representation Eq. ( p5| ) is completely 
general if we interpret T as a function of oj which in linear- 
response theory is related to the neutron spin's "memory 
function" 0]. In our heuristic description, however, we 
will use a constant value for T which is fixed by the nor- 
malization requirement Eq. (|l0|). 

In order to calculate r CT and s(x) in a dilute neutron 
medium we model the nn interaction by one-pion ex- 
change in Born approximation, an approach which has 
been common practice for SN and neutron-star physics 
since Friman and Maxwell's seminal paper |l3| and which 
is further justified in Ref. fji| ]. Further, we take the neu- 
trons to be nondegenerate which is not a bad approxi- 
mation during the early phases of SN core cooling. Fi- 
nally, we neglect the mass in the pion propagator which 



is also a reasonable approximation for the large momen- 
tum transfers in typical nn collisions in a SN core. All of 
these approximations go in the same direction of some- 
what overestimating the nn spin interaction rate. We 
also ignore static spin-spin correlations which could, in 
principle, both enhance or diminish our results. 

Within this framework the spin-fluctuation rate is ex- 
plicitly found to be 



T a =A^aln n T 1 ' 2 r 



-5/2 



N 



(26) 



where a n = (f2m,N /m^) 2 /An w 15 with / w 1 is the 
pion-nucleon "fine-structure constant," n n is the neutron 
density, and m^v the nucleon mass. Numerically we find 



where pi± e 
over, one finds 111 




i.6p 14 T 10 1/2 



cm and T- 



10 



(27) 

T/lOMeV. More- 



s(x) 



dv e 



<(o, 



2(2v- 



\J v(v + x) 
- log 



\fv + X + 



\/v + x — 



(28) 



an expression which indeed fulfills the detailed balance 
requirement Eq. (E4f) and which is smooth at x = with 
the derivative s'(0) = 1/2 (Fig. |). We will use a simple 
analytic approximation to this integral 1141] 



s{x) 




12+3 

7T 



-1/12 N 



-1/2 



(29) 



which reproduces the correct limiting behavior for i»l 
and for x — where it also has the correct derivative. It 
deviates from the true value by no more than 2.5% any- 
where. For x < we use s(x) — s(—x) e x in accordance 
with detailed balance. 




10.0 



FIG. 1. The function s(x) as denned in Eq. (|28|). The 
analytic approximation Eq. (E9() is identical to within plotting 
accuracy. 
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10.0 15.0 20.0 25.0 30.0 

Y,=r„/T 

FIG. 2. r appearing in Eq. to normalize Scr(uj) as a 
function of r CT . 




10.0 15.0 20.0 25.0 30.0 

Y„ = r „ / T 

FIG. 3. /-sum of S<t(lj) as defined in Eq. ( |30| ) as a function 

of r CT . 




10.0 15.0 20.0 25.0 30.0 

y o = r o / t 

FIG. 4. Static long- wavel eng th limit of the magnetic sus- 
ceptibility according to Eq. (h9h as a function of r CT . 



Our semiheuristic toy model is thus completely de- 
fined. In Fig. ^ we show the "downstairs T" of Eq. ( p5|) as 
a function of T a such that the normalization requirement 
Eq. (|l(]) is obeyed. By construction we have V = T a for 
— > with smaller values for a larger T a . This reduc- 
tion is mostly due to the detailed-balance behavior which 
suppresses the classical structure function for negative u>. 

Further we consider the /-sum which is for our present 
model 



+ 30 



(kt) 
2^ 



UlS a {uj) = T 



tlx 



2tt x 2 + 7 2 /4 



s(x), (30) 



where 7 = T/T. As claimed before it is equal to the 
spin fluctuation rate times a factor of order unity which 
is shown in Fig. | as a function of T a . 

Finally we show in Fig. |] the static long-wavelength 
susceptibility in units of the Pauli susceptibility for our 
toy model according to Eq. ([H)]) . It is a slowly decreasing 
function of the spin fluctuation rate. 

The overall spin-density structure function S a (uj) in 
our toy model is shown in Fig. || for several values of 
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FIG. 5. Dynamical spin-density structure function S a (u>) 
in our heuristic model Eq. (|2|) for 70- = 3, 10, and 30. 

Next we study the photon dispersion relation implied 
by our model. In Fig. ^ we show the long-wavelength 
limit x'(w) = ^refr — ^ m un hs of the Pauli susceptibility 
as a function of x = lu/T for several values of 7 CT . The 
overall behavior i s exa ctly as expected from our general 



discussion in Sec. IIP. To see the large-w behavior more 



clearly we show in Fig. [?] the equivalent quantity 



in units of XPauiiT a T. As predicted, 



'err 



'err 



approaches an 



asymptotic value which is independent of frequency and 
which is of order XPauirr CT T. Of course, for small uu the 
squared "effective mass" m 2 s begins at negative values. 
However, for all frequencies and all values of 7^ we find 
< m 2 where the latter is the asymptotic value 



that \rn 
for uj 
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FIG. 6. Dispersive part \ °f the magnetic susceptibility 
in our toy model in units of XPauii- 
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FIG. 7. Effective photon mass in our toy model as a func- 
tion of photon frequency. 



In order to arrive at a numerical estimate for the 
magnetically induced photon transverse mass we write 
the nucleon magnetic moments in the usual form //jv = 
KAre/2mjv with to^v the nucleon mass and K n = —1.91 
and k p — 2.79. For simplicity we treat all nucleons as if 
they were neutrons which implies 



XPauii = 4.4x 10- 3 pu/T] 



10. 



(31) 



where p u = p/10 14 gem" 3 and Tin = T/lOMeV. 
With the perturbative estimate Eq. (P7[) for the spin- 
fluctuation rate we then find 



in 



,r y I magnetic 



(XPauhTT,) 1 / 2 « 1.9 MeV puT^*. (32) 



In a SN core we have densities of up to 10 15 g cm -3 and 
temperatures of up to 30-60 MeV, implying a perturba- 
tive spin-fluctuation rate far in excess of the temperature, 
i.e. 7o- = 10-100. It has been argued that in a SN core the 
true r<j cannot exceed a few times T ||,|l^|- Therefore, 
we have probably overestimated the transverse photon 
mass by at least a factor of a few. 



We next consider the corresponding quantity caused by 
the interaction with electrons. One finds jl6| m 2 = \lo p 
with ojp the plasma frequency. For relativistic degenerate 
electrons it is top = (4a/3iv)pp e with pp :C the electron 
Fermi momentum so that 



7 I electrons 



2a 2 
— Pf, ( 

7T 



Numerically, this corresponds to 



7 I electrons 



= 16.3 MeV K 



1/3 1/3 
Pl4 



(33) 



(34) 



Evi- 



where Y e is the number of electrons per baryon. 
dently in the center of a SN core with pi4 ps 8, Tin ps 4 
and initially Y e ps 0.3 the magnetic moment contribution 
could be almost as large as the electronic term. However, 
because we have probably overestimated the magnetic 
term by a factor of a few the electrons still dominate. 



IV. DISCUSSION AND SUMMARY 

We have calculated the photon refractive index due to 
the interaction with the magnetic moments of the nu- 
cleons. For simplicity we have limited our discussion 
to nondegenerate neutrons. In the collisionlcss limit the 
forward-scattering amplitude vanishes identically so that 
the neutron magnetic moments alone do not cause any 
deviation of the photon dispersion relation from the vac- 
uum behavior |3j. However, because of strong neutron 
spin interactions the collisonless limit is far from justified 
in a SN core. On the basis of the fluctuation-dissipation 
theorem and the Kramers-Kronig relation we have de- 
rived a general expression for the photon refractive index 
in terms of the dynamical neutron spin-density structure 
function S <j {uj 1 k). In an interacting medium it is a broad 
function of u, in contrast to the collisionless limit where 
it is proportional to 8(lo). 

We have found that for w « r a (the neutron spin 
fluctuation rate) the "effective photon mass" m^ s begins 
with negative values — XPauii^ 2 m terms of the Pauli sus- 



ceptibility of the neutron ensemble. However, as shown in 
Fig. |?] this function quickly turns around and then grows 
asymptotically to a positive value m? ( m XPauii2T CT . In 
absolute terms this "transverse photon mass" is much 
larger than the maximum excursion of rn^ s to negative 
values. 

A numerical comparison for conditions relevant for a 
SN core reveals that the transverse photon mass caused 
by the neutron magnetic moment tends to be much 
smaller than that caused by the electron plasma effect, 
except for extreme densities and low electron fractions 
where the magnetic term may actually compete with the 
electronic one. A numerically accurate comparison is 
not possible because the neutron dynamical spin-density 
structure function is not known in any detail. We have 
only performed a relatively schematic estimate which in- 
volved many simplifying assumptions. However, it still 
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appears safe to conclude that the negative magnetic m^g 
at small frequencies cannot compete with the electronic 
plasma effect. This indicates that the total m^ ff is al- 
ways positive, i.e. the photon refractive index is always 
less than 1 and it is reasonably well estimated by the 
electronic plasma effect. This implies that the Cherenkov 
processes v <-* remain forbidden. 

A more quantitative analysis than has been presented 
here requires a better understanding of the dynamical nu- 
cleon spin-density structure function, or more precisely, 
of the dynamical spin and isospin susceptibilities of a 
hot and dense nuclear medium. We stress that the uj 
dependence is crucial for the photon dispersion relation 
as well as the neutrino opacities j|, the static suscepti- 
bilities alone which have sometimes been studied in the 
literature are not enough. 
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